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ABSTRACT 

We present resolved images of the dust continuum emission from the debris 
disk around the young (80-200 Myr) solar-type star HD 107146 with CARMA at 
A1.3mm and the CSO at A350 /im. Both images show that the dust emission 
extends over an ~10" diameter region. The high resolution (3") CARMA image 
further reveals that the dust is distributed in a partial ring with significant de- 
crease in flux inward of 97 AU. Two prominent emission peaks appear within the 
ring separated by ~140° in position angle. The morphology of the dust emission 
is suggestive of dust captured into a mean motion resonance, which would imply 
the presence of a planet at an orbital radius of ~45-75 AU. 

Subject headings: circumstellar matter-planetary systems-stars: individual (HD 
107146) 

1. Introduction 

Photometric surveys with IRAS, ISO, and the Spitzer Space Telescope have identified 
hundreds of main-sequence stars where the observed flux densities at mid-infrared wave- 
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lengths are brighter than the stellar photosphere (for a review, see |Meyer et aL|2007l ). This 
"excess" infrared emission is generally attributed to small dust grains that are produced 
when planetesimals collide and fragment into smaller debris. 

The gravitational interaction of any planets with the debris dust produces gaps and 
peaks in the dust spatial distribution that provides clues to the underlying planetary system 
( Liou fc Zook| l~999 ) . However, the dust emission in most debris systems remains unresolved, 
and the dust spatial distribution is often inferred by modeling the composite spectral en- 
ergy distribution. In practice such analyses yield ambiguous results, since simplified disk 
geometries are assumed, and the inferred dust spatial distribution depends sensitively on the 
assumed grain composition, porosity, and size. 

A few nearby debris disks have been resolved, providing direct information on the dust 
spatial distribution. Scattered light images at visible and near-infrared images often show 



smooth azimuthal distributions, frequently distributed in ring-like geometries (Kalas, Liu, & 



Matthews 2004; Kalas, Graham, & Clampin 2005 Schneider et al. 2006). Resolved millimeter 



wavelength images also show rings of material, but with more clumpy distributions of debris 



emission (Holland et al. 1998 Greaves et al. 1998 2005 Koerner, Sargent, & Ostroff 2001 



Wilner et al. 2002; Maness et al. 2008). These clumps have been interpreted as signposts 



of planetary systems where the gravitational interactions trap the dust in mean motion 



resonances (Wilner et al. 2002 Wyatt 2003) 



The debris disk around the G2V star HD 107146 (Jaschek 1978) has received special 



interest since it is a rare example of a resolved debris disk around a young (80-200 Myr; 



Williams et al. 2004) solar analog. At a distance of 28.5 pc (Perryman et al. 1997), the 



disk has been resolved in near-infrared scattered light (Ardila et al. 2004) and in the mil 



limeter thermal continuum (Williams et al. 2004; Carpenter et al. 2005). We present new 



(sub-) millimeter- wavelength images of the HD 107146 debris disk obtained with CARMA 
(Combined Array for Research in Millimeter-wave Astronomy) and the CSO (Caltech Sub- 
millimeter Observatory). The debris disk is clearly resolved in both datasets, with the 
CARMA image revealing a complex morphology not seen previously. We describe the obser- 
vations in £[2j present the resolved images in £|3j and discuss the implications of these results 
in m 



2. Observations 



Continuum images of HD 107146 were obtained at wavelengths of 1.3 mm and 350 /im 
using CARMA and the CSO, respectively. We adopted equatorial coordinates of (a, 5) = 
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(12:19:06. 50, +16:32:53. 87) for equinox and epoch J2000, with a proper motion of (fi a ,f^s) 



;-175.65, -148.28) mas yr" 1 (Perryman et al. 1997). 



CSO observations were carried out between UT April 17 and 21, 2005, using the Sub- 



millimeter High Angular Resolution Camera II (SHARC II; Dowell et al. 2003) with the 
350 fim filter. The median zenith opacity at 225 GHz was 0.04, which equates to an opacity 
of ~1 at 350 /im. Telescope pointing was checked on Callisto every 30 minutes, which was 
20° away from HD 107146 at the time of the observations. Pointing was also checked on 
the fainter source 3C273, which is 15° from HD 107146, every hour. The pointing offsets 
as a function of time and elevation were interpolated in offline data processing and applied 
before co-adding the HD 107146 data. Absolute flux calibration was set from observations of 
Neptune assuming a total flux density of 92.6 Jy. We estimate a la calibration uncertainty 
of 15% from repeated observations of calibrators over the entire observing run. Images 



were produced using the Comprehensive Reduction Utility for SHARC II (CRUSH; Kovacs 



2008). The "deep" data reduction mode in CRUSH was applied to both the calibrators and 



HD 107146. 

CARMA observations were conducted on UT 2007 August 21, 2008 September 12, and 
2008 September 14 in E-configuration, and on four days between UT March 3 and March 
7, 2008 in D-configuration. The E- and D- configurations span baselines of 6-66 m and 
1 1-148 m, respectively. Observations were conducted with the local oscillator (LO) set to a 
frequency of 227.25 GHz. Three correlator bands of 468 MHz bandwidth each were placed 
at intermediate frequencies (IF) of 2.25, 2.75, and 3.25 GHz to provide a total continuum 
bandwidth of 2.8 GHz after combining the upper and lower sidebands. The quasar 3C273 was 
observed approximately every 20 minutes for amplitude, phase, and passband calibration. 
For the March and September 2008 observations, pointing was updated every ~30 minutes 
using optical reference pointing (Corder, Carpenter, & Wright 2008, in preparation). 



All reduction and imaging were performed using MIRIAD (Sault, Teuben, & Wright 



1995). Data were flux calibrated by observing Mars or MWC349. For MWC349, we as- 



sumed a flux density of 1.69 Jy following the calibration adopted at the PdBlj^] with an 



assumed uncertainty of 10% (Altenhoff, Thum, & Wendker 1994). The night-to-night RMS 



repeatability of the measure flux density of 3C273 is 2.4%. We adopt a net calibration 
uncertainty of 15%. The image was formed with natural weighting of the uv visibilities 
and a Gaussian taper of FWHM=1.6", and then deconvolved with the point spread function 
(i.e. the "dirty" beam) using a hybrid Hogbom/ Clark/ Steer algorithm ( Hogbom]|1974 ; Clark 



1980 Steer, Dewdney, & Ito 1984). Since the emission is contained to within a radius of 



^ttp:// ww.iram.fr/IRAMFR/IS 
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<5" from the phase center, no primary beam corrections were applied. The final image has 
a resolution of 3.2"x2.7" and a RMS noise of 0.35 mJy beam -1 . 

3. Resolved Images 

Figure [T] presents the CSO 350 /an continuum image of HD 107146. The RMS noise in 
the image is 6mJy beam -1 , and the peak flux density (181 mJy beam -1 ) is detected with a 
signal to noise ratio of 30. The integrated 350 /im flux density obtained by fitting an elliptical 
Gaussian to the image is 319±45mJy, where the la uncertainty includes the statistical 
(6mJy) and calibration (45mJy) uncertainties added in quadrature. The centroid of the 
350/im continuum emission is offset from the stellar position by (Aa, AS) = (0.5",-1.0"); 
this offset is consistent within the pointing accuracy of the CSO telescope. No significant 
azimuthal structure in the 350 /im emission is detected. 

Figure [2] shows the azimuthally-averaged radial profile of the 350 /im emission for HD 107146 
(solid curve) and Callisto (dashed curve). Callisto had an angular diameter of 1.49" at the 
time of the observations and is considered a point source for this analysis. The measured 
FWHM of the emission profile on Callisto is 9.8"±0.03". By contrast, the FWHM size of 
the HD 107146 debris disk measured from an elliptical Gaussian fit to the 350 /im image is 
13.5" x 12.6". The deconvolved image size of HD 107146 using the observed Callisto image 
to represent the beam profile is (8.9" ±0.6") x (8.2" ±0.5"). The 350 /tm continuum emission 
toward HD 107146 is clearly resolved. 

Figure [3] presents the CARMA 1.3 mm continuum image of HD 107146 at 3" resolution. 
The 1.3 mm continuum emission extends over a ~10" region. The integrated flux density 
obtained by summing the emission within a 12" x 12" box centered on the stellar position is 
10.4 ± 1.4 mJy, with an additional 15% calibration uncertainty (§]2]). 

The 1.3 mm emission is resolved into a partial ring that encircles the stellar position. 
The locations of the emission peaks are robust to various calibration schemes and data 
reduction approaches, but the magnitude of the peaks and the gap in the northwestern part 
of the ring are somewhat variable with such changes. The two brightest positions in the ring 
are located nearly equidistant from the stellar position on opposite sides of the star. The 
northeast peak is located 3.3" from the star at a position angle (east of north) of PA=48°, 
and the southwest peak is offset by 3.5" from the stellar position at PA=190°. The projected 
distance of the dust peaks from the star is ~97AU. 



The clumpy structure observed in the CARMA 1.3 mm continuum image is in contrast 
to the smooth, continuous ring observed in scattered light at optical wavelengths (Ardila 
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Fig. 1. — CSO 350 /iin continuum image of the HD 107146 debris disk. Contours begin at 
3a with increments of 5a, where a = 6mJy beam -1 . The filled star at (0,0) indicates the 
stellar position. The filled circle in the lower left corner represents the 9.8" resolution of the 
observations measured on the Callisto image. 
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Fig. 2. — Radial profile of the CSO 350 /zm continuum emission for the HD 107146 debris 
disk (solid curve) and Callisto (dashed curve). Both radial profiles were measured on coadded 
images of the sources. 



et al. 2004). Further, the scattered-light image peaks at a larger orbital radius (130 AU) 



than the dust continuum image. The differences in the radial distribution is qualitatively 
consistent with the notion that small grains, which are efficient at scattering, are ejected 
by radiation pressure from the system, while the larger grains are found at smaller radii 



(Takeuchi & Artymowicz 2001 Wyatt 2006). 



The morphological differences of the 350 /im and 1.3 mm images can be attributed to 
resolution in that if we produce the 1.3 mm CARMA image with a 8" FWHM taper to yield 
a 9" synthesized beam, the resulting image is centrally peaked on the stellar position and 
is similar in size to the 350 /im image. In detail, other published, resolved, (sub-)millimeter 



images are inconsistent with the CSO and CARMA images. Williams et al. (2004) report 



that the 450 /im emission at 8" resolution is offset from the stellar position by 4.4". Since the 
350 /im (see Fig. Q and the 850 /zm image ( |Williams et al.|2004[ ) are centered on the stellar 
position to within 2", we assume the offset in the 450 /im image is a pointing error. A more 
significant discrepancy is seen between the CARMA 1.3 mm image and the OVRO 3 mm 
image ( |Carpenter et~aT]|2l)05~| ). The 3 mm emission was resolved at 4.5'" x 4.0" resolution 
with a FWHM size of (6.5"± 1.4) x (4.2" ± 1.3"). We have smoothed the 1.3 mm image to the 
resolution of the 3 mm image, and the ring structure is not apparent at the coarser resolution. 
However, the apparent ellipticity in the 3 mm image is nearly orthogonal to the clumps in 
the 1.3 mm image. Since the 1.3 mm image has higher signal to noise, finer resolution, and 
more baselines for better imaging, we view that image as a more reliable tracer of the debris 
structure. 



4. Implications 



Assuming gas drag is negligible, the orbital lifetime of dust grains in an optically thin 
disk is limited by Poynting-Robertson drag, collision grinding of particles combined with 



radiation blowout of small grains, and stellar- wind drag. Dominik & Decin (2003) have 



shown that for debris disks detectable with current instrumentation, collisional grinding of 
particles to the radiation blowout size dominates over Poynting-Robertson drag in removing 



mass from the debris system. Stellar- wind drag may also be important (Gustafson 1994), 



but the magnitude of this effect is poorly constrained at young ages (Jura 2004; Wood et 



al. 2005). 



If a collisional cascade has produced the dust particles in the HD 107146 debris disk, 
the expected particle size distribution is n(a) oc a~ 3 5 (Dohnanyi 1969) for grain sizes larger 
than the radiation blowout size (radius ~0.5/aul for a solar-type star). In practice, the 
truncated distribution will produce a "wavy" pattern of particle sizes superimposed on the 
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Fig. 3.— CARMA 1.3 mm continuum image of the HD 107146 debris disk. The FWHM of 
the synthesized beam (3.2" x 2.7") is shown in the lower left corner. Solid contours begin at 
2a with increments of la, where a = 0.35 mJy beam -1 . Dashed contours begin at —2a; no 
< — 3a pixel values are present within the region shown. The filled star at (0,0) indicates the 
stellar position. The "dirty" beam has been deconvolved, but no primary beam corrections 
have been applied. 
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power-law (e.g. Thebault & Augereau 2007), but we neglect that complication. For silicate 



grains ( Weingartner & Draine 2001) larger than the blowout size, half of the flux density 
at a wavelength of 1.3 mm will be emitted by grains larger than ~0.5mm in radius for a 
q-3.5 g j ze distribution. The ratio of radiation to gravitational forces {(3) for these particles 
is /5~0.001, indicating radiation pressure has negligible effect on their dynamics. 

By the above arguments, the millimeter-wave emission traces predominantly millimeter- 
sized particles where the dynamics are controlled by particle collisions and gravity. This is 
consistent with clumpy structure observed in the debris disk, as radiation forces will tend 
to produce a smooth spatial distribution of particles. The millimeter-wave emission should 
then trace the location of the planetesimals in the HD 107146 disk ( Wyatt||2006 ). The radius 
of the observed debris emission implies the presence of planetesimals at an orbital radius 
of ~97AU. Moreover, the dual emission peaks in the CARMA image presents intriguing 
possibilities for the presumed planetary system around HD 107146. The dust emission peaks 
are distributed nearly (but not precisely) symmetrically around the star, and may be caused 
by dust trapped in mean-motion resonance with an orbiting star or planet. In fact, the image 



of HD 107146 looks quite similar to the /3~0.002, 3:2 resonance image of Wyatt (2006). In 



detail, the two main dust peaks and the star are not co-linear, which may imply a planet of 
moderate eccentricity (e ~0.5) has trapped the dust in a 3:1 resonance (Kuchner & Holman 



2003). Given that the clumps are located at 97 AU, such resonances would imply a planet 
between ~45 and 75 AU. 



Metchev & Hillenbrand (2003) obtained adaptive optics images at A2.2/zm to search for 



companions to HD 107146. They place an upper mass limit of 10 Mj to any companions at 
74 AU. Apai et al. (2008) report that a similar limit of 10 Mj persists inward to ~25AU. 



Thus these observations rule out the presence of a stellar companion at such a large radius, 
and suggest the structure in the debris disk is produced from the gravitational affects of a 
planetary system. 

Could a planet have formed at ~45-75AU that produced resonant structure in the 



debris disk? Kenyon & Bromley (2004a) have simulated the collision growth of planets and 
subsequent production of debris and find that planets on the size of ~2000km (~0.01M e ) 
are needed to excite the collisional cascade. The timescale to produce bodies of this size at 
45-75 AU from collisional growth of planetesimals is ~50-225 Myr ( Kenyon fc Bromley|2004a ) 
for planetesimal with low eccentricity (e < 0.02) and a surface density profile appropriate 
for the minimum mass solar nebula ( Weidenschilling 1977). For moderate eccentricities 
(e ~0. 02-0. 04), the timescales are 2 to 4 times longer. Similarly, Kenyon & Bromley (2004b) 



simulated the collisional growth of objects at 40-47 AU and found that 1000 km sized objects 
can form in 10-50 Myr, and double in mass every 100 Myr to 1 Gyr. Thus the timescale to 
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form ~2000km sized bodies, which are needed to initiate the collisional cascade, is roughly 
consistent with the age estimates of HD 107146. 

However, the planet is likely much larger than 2000 km in radius if the dust is trapped 
in a resonance. Roques et al.l (1994) conducted numerical simulations to show that planets 



more massive than ~5M e (~ll,000km in radius assuming a volume density comparable to 
the Earth) are needed to trap particles in outer mean motion resonances. The timescale 



to grow such objects by collisional growth is more than 1 Gyr at radii >45AU (Kenyon & 



Bromley 2004b), which is substantially older than HD 107146. 



The presence of a massive planet at large radii then requires formation mechanisms 
beyond collisional aggregation of planetesimals. For disk-to-star mass ratios of >0.1 in the 
T Tauri phase, marginally consistent with observations ( Andrews fc Williams|[2007 ), planets 
may form by gravitational instability (e.g. Boss|200"8 and references therein). Planets formed 



by this mechanism tend to be quite massive (lMj or larger). Rafikov (2005) argues that 
the cooling timescales are inconsistent with formation via this method unless the planet is 
both massive, ~5-10Mj, and distant, ~100AU, which is consistent with the distance and 
upper mass limit described above. Another possibility is that the planet formed in the inner 
regions and migrated outwards, which may trap planetesimals and dust in resonant orbits 



(Wyatt 2003). Such a mechanism has been used to explain the structure and dynamics of 



the Kuiper Belt ( Malhotra|1993 , 1995). Detailed modeling is needed to gain more insights on 
the planetary system around HD 107146 and if migration can explain the clumpy structure 
observed in the surrounding debris disk. 

In summary, the young solar analog HD 107146 is surrounded by a clumpy debris disk 
similar to those seen around Vega and eEridani. The structure is suggestive of dynamical 
influence from a planet. The likely size and location of the presumed planet are inconsistent 
with formation by purely collisional aggregation of planetesimals and appeals to migration 
or formation via gravitational instability must be made. 
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